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Phagocytosis is the fundamental ceUular process by which eukaryotic ceUs bind and engulf particles 
by their cell membrane. Particle engulfment involves particle recognition by cell-surface receptors, 
signaling and remodeling of the actin cytoskeleton to guide the membrane around the particle 
in a zipper-like fashion. Despite the signaling complexity, phagocytosis also depends strongly on 
biophysical parameters, such as particle shape, and the need for actin-driven force generation remains 
poorly understood. Here, we propose a novel, three-dimensional and stochastic biophysical model 
of phagocytosis, and study the engulfment of particles of various sizes and shapes, including spiral 
and rod-shaped particles reminiscent of bacteria. Highly curved shapes are not taken up, in line 
with recent experimental results. Furthermore, we surprisingly find that even without actin-driven 
force generation, engulfment proceeds in a large regime of parameter values, albeit more slowly and 
with highly variable phagocytic cups. We experimentally confirm these predictions using fibroblasts, 
transfected with immunoreceptor Fc7RHa for engulfment of immunoglobulin G-opsonized particles. 
Specifically, we compare the wild-type receptor with a mutant receptor, unable to signal to the 
actin cytoskeleton. Based on the reconstruction of phagocytic cups from imaging data, we indeed 
show that cells are able to engulf small particles even without support from biological actin-driven 
processes. This suggests that biochemical pathways render the evolutionary ancient process of 
phagocytic highly robust, allowing cells to engulf even very large particles. The particle-shape 
dependence of phagocytosis makes a systematic investigation of host-pathogen interactions and an 
efficient design of a vehicle for drug delivery possible. 



Background 



Phagocytosis is the ancient, evolutionarily conserved 
process by which eukaryotic cells bind, engulf, and de- 
stroy particles and cells larger than 0.5/im in diameter 
[1-3]. The importance of phagocytosis is derived from its 
two main functions: (1) a feeding mechanism in single- 
cell organisms [4], and (2) the clearance of pathogens, 
apoptotic and senescent cells from our body by immune 
cells [5, 6]. As part of our immune defense, phagocytosis 
is mainly performed by professional phagocytes, includ- 
ing macrophages, neutrophils, and dendritic cells. Initi- 
ation of phagocytosis occurs with recognition of the tar- 
get particle either directly or via an opsonising molecule. 
For instance the Fc portion of immunoglobulin G (IgG) 
is recognized by the cell-surface receptor Fc7RIIa [7, 8]. 
Ligand-receptor binding triggers intracellular signaling 
[3, 8, 9], resulting in remodeling of the actin cytoskeleton 
[10, 11] and coherent growth of cell membrane around the 
particle to form the phagocytic cup [1, 2]. Eventually, the 
leading edge of the growing cup closes, and a membrane 
vesicle enclosing the particle (phagosome) moves inside 
the cell. Subsequently, the phagosome fuses with vesicles 
containing enzymes [12, 13], acids [14], and oxygen radi- 
cals [15, 16] to destroy the particle. 

The biochemical pathways involved in phagocytosis are 
complex. Dozens of cell-surface receptors contribute to 



the recognition of a large variety of ligand molecules and 
subsequent particle engulfment [3, 8, 17]. The Fc7 re- 
ceptor (FC7R) [18] and complement receptor 3 (CR3) 
of the integrin receptor family [19] are the most widely 
studied and understood receptors involved in phagocyto- 
sis. Fc7R-mediated phagocytosis proceeds through mem- 
brane protrusions and leads to thin cups [20, 21], whereas 
in CR3-mediated phagocytosis, particles appear to sink 
into the cell [22, 23]. Spreading of the cell membrane 
over the particle involves actin-driven cell-shape changes 
similar to the processes involved in cell migration and 
adhesion [7, 24, 25, 27]. Specifically for FC7R, binding 
to an IgG-opsonized particle results in receptor cluster- 
ing and recruitment of small GTPases of the Rho fam- 
ily, which, via proteins of the WASP family, activate the 
Arp2/3 complex [1, 3]. The latter promotes branching of 
actin filaments, leading to an increase in the number of 
uncapped ends and to an isotropic growth of the actin 
network [8, 28]. Additionally, the phagocytic cup has 
been shown to be enriched in gelsolin [29-31], coronin 
[4], and other regulators of actin polymerization. All in 
all, this complex signaling pathway involves 100-1000 dif- 
ferent types of molecules [3, 32], rendering mathematical 
modeling at the molecular level impossible. 

Despite the huge biochemical complexity, the engulfment 
process shows a strong dependence on simple biophysi- 
cal parameters. First, it relies on the availability of ex- 
tra membrane at the phagocytic cup [33, 34], provided 
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by delivery of membrane vesicles [35] or unwrinkling of 
membrane folds [36, 37]. Second, completion of phago- 
cytic uptake depends on the shape of the particle and, 
interestingly, on the initial orientation of the particle 
on the cell surface [20, 38]. For instance, experiments 
demonstrate that elongated spheroid polystyrene parti- 
cles coated with IgG are more efficiently engulfed when 
presented to the phagocyte with their tip first. Third, a 
recent study by one of the authors demonstrates that the 
biophysical requirements for phagocytosis lead to either 
complete phagocytosis or stalled cups due to the presence 
of a mechanical bottleneck [39]. Interestingly, the same 
study shows that engulfment appears to even proceed in 
cells treated with (modest amounts of) cytochalasin D, 
an inhibitor of actin polymerization, indicating that bio- 
chemical pathways may not always be necessary for this 
initial stage of phagocytosis. 

The mechanism of phagocytosis is only partially un- 
derstood, with key insights provided more than three 
decades ago. In the 1970's, Griffin and his collabora- 
tors [40, 41] demonstrated that incomplete coating of 
particles with ligand results in only partial uptake. This 
indicated that phagocytic uptake occurs via successive 
zipper-like ligand-receptor binding (Figure lA), and not 
by an all-or-nothing mechanism triggered at the onset 
of phagocytosis. The zipper mechanism is the underly- 
ing assumption in a number of recent modeling works 
in phagocytosis [39, 42] and endocytosis [43-45], mainly 
addressing the infiuence of the cell-membrane tension 
and ligand-receptor bond density on engulfment. De- 
spite the general acceptance of the zipper mechanism, 
many of its biophysical requirements are insufficiently 
understood. Questions, so far unanswered, include what 
the energetic requirements of the zipper mechanism are, 
specifically what role actin polymerization plays in its 
progression during phagocytosis, and also whether the 
zipper mechanism can explain the particle-shape depen- 
dence of phagocytosis. Previous models were unable to 
fully address the particle-shape dependence, as they as- 
sume rotational symmetry around the axis connecting 
cell and particle. Additionally, large particle-to-particle 
variation in cup growth [39] and cell-to-cell variation in 
the related process of endocytosis [46] point towards the 
importance of stochasticity during the uptake, not cap- 
tured in previous deterministic approaches. 

Recent experiments provide new insights into the bio- 
physical mechanism for driving the membrane around 
the particle, suggesting a ratchet-type mechanism. Once 
started, phagocytosis progresses unidirectionally and ir- 
reversibly [37]. This irreversible membrane progression 
is further supported by the loss of lipid and protein mo- 
bility at the phagocytic cup, observed using fiuorescence 
recovery after photobleaching (FRAP) [47]. While sev- 
eral models proposed mechanisms of force generation by 
actin polymerization (see [48] and references therein), re- 
cent experiments based on fiuorescent speckle microscopy 



demonstrate that actin does not directly push the mem- 
brane outwards. Instead, by filling gaps provided by 
membrane fiuctuations (or other types of membrane 
movement), actin polymerization prevents the membrane 
from moving backwards like a ratchet [42]. The relevance 
of such Brownian ratchets in biology has previously been 
emphasized [49, 50]. The question is if a ratchet mech- 
anism, together with energetic restrictions in membrane 
bending and stretching, can naturally lead to phagocytic 
uptake and account for the shape-dependence of phago- 
cytosis. 

In this work, we propose a ratchet-like biophysical model 
for the zipper mechanism. This model differs from pre- 
vious works in that it is, to our knowledge, the first fully 
three-dimensional stochastic model of phagocytic engulf- 
ment. Specifically, thermal membrane fiuctuations, as- 
sumed to play a major role in our model, provide the en- 
ergy source to locally deform the membrane and to build 
further ligand-receptor bonds for zippering the mem- 
brane around the particle. Actin polymerization makes 
ligand-receptor bonds effectively irreversible, i.e. rein- 
forced and stabilized for a significant amount of time. 
To investigate the role of actin, we compare cup progres- 
sion for the regular active zipper with a passive zipper 
model in which ligand-receptor binding remains specific 
and strong but reversible due to the absence of actin poly- 
merization. 

Interestingly, we find that the passive zipper also leads 
to engulfment of small particles, rendering phagocyto- 
sis highly robust. However, such passive engulfment is 
generally slower and produces much more variable phago- 
cytic cups than the active zipper. Furthermore, our com- 
puter simulations lead to successful phagocytic engulf- 
ment in a broad range of parameters values, including 
different particle sizes. For non-spherical particles, com- 
pletion of engulfment depends strongly on particle shape 
and orientation. Our model further predicts that cup 
shape invariably depends on membrane biophysical pa- 
rameters, in particular surface tension and cell- volume 
constraint. 

To test the predicted difference between the active and 
passive zippers, we experimentally implement the two dif- 
ferent types of zippers using COS-7 fibroblasts which, 
after transfection with GFP-tagged Fc7 receptor, phago- 
cytoze IgG-coated polystyrene particles. Specifically, we 
performed phagocytic assays under three different condi- 
tions: (1) cells expressing wild- type FC7R for the active 
zipper (WT-FC7R), (2) cells expressing a signaling-dead 
mutant receptor (Y282F/Y298F-Fc7R), which specifi- 
cally binds IgG ligand but is unable to signal to the 
actin cytoskeleton [39, 51, 52], and (3) cells express- 
ing WT-FC7R and treated with cytochalasin D (WT- 
Fc7R+CytoD). The last two conditions represent two 
versions of the passive zipper due to the absence of actin 
polymerization in phagocytic cups. To compare with our 
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FIG. 1: (A) Schematic of the zipper mechanism. Cell-membrane receptors (green dots) and ligands on the particle surface 
(blue dots) are sequentially engaged in bond formation, resulting in progression of engulfment with time (from left to right). 
(B) Schematic of our ratchet model, {left) A random membrane fluctuation (blue) far from the particle is unable to trigger 
ligand-receptor binding and signaling. Therefore it is not supported by remodeling of the actin cortex (straight red line), and 
the membrane may move backwards at a later time, {right) A membrane fluctuation near the particle leads to ligand-receptor 
binding, resulting in signaling and actin polymerization (red diamonds) . Consequently, the actin cortex is deformed to support 
the membrane fluctuation, which makes the membrane move irreversibly for zipper progression. Energetic costs and gains of 
membrane fluctuations used in model are shown as well. 



model, we systematically analyze confocal microscopy 
images, and quantitatively estimate cup variability for 
the three different conditions using small (1.5 /im radius) 
and large (3 /im radius) particles. Consistently with our 
simulations, phagocytic cups develop more slowly and are 
significantly more variable in the absence of actin poly- 
merization. Our results provide new insights into the ro- 
bustness of phagocytosis, as well as the role of bacterial 
cell shape in host-pathogen interactions. 



Results and Discussion 

Ratchet model for the zipper mechanism 

Our model is based on the following experimental ob- 
servations. Engulfment of quasi-spherical particles by 
neutrophils progresses continuously without significant 
pause or reversal, indicating that ligand-receptor bind- 
ing is essentially irreversible [37]. This irreversibility is 
further supported by FRAP and single-molecule exper- 



iments, which show that lipids and proteins in phago- 
cytic cups, as well as ligand-bound Fc-receptors are im- 
mobilized in an actin-dependent fashion [47, 53]. Addi- 
tional support for the notion of irreversible uptake was 
recently determined in a related context [54] . Fluorescent 
speckle microscopy of actin flow and image analysis dur- 
ing cell migration show that the membrane at the leading 
edge protrudes first, followed by actin polymerization to 
fill the gap between the membrane and the actin cor- 
tex. Such actin polymerization is mainly restricted to 
the leading edge due to signaling by receptors and/or lo- 
calization of small GTPases of the Rho family [8]. The 
role of actin polymerization in phagocytosis is hence to 
stabilize ligand-receptor bonds and to rectify membrane 
movements in a ratchet-like fashion, leading to unidirec- 
tional movement of the leading edge of the engulfing cell. 

Figure lA introduces the general concept of the zipper 
mechanism in phagocytosis, and Figure IB summarizes 
our ratchet model for this mechanism. The cell mem- 
brane and the actin cortex are described by a Helfrich- 
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FIG. 2: (A) Influence of surface tension and volume constraint. Simulated phagocytic cup for {left) a large volume constraint 
(kp = 10~^ pN/im~^) and a relatively low surface tension (a = 6.2 x 10~^ mNm~^), and (right) a high surface tension 
(cr = 6.2 X 10~^ mNm~^) and a nearly unconstrained cell volume (kp = 3.8 x 10"'' pN/im~^). Both sets of parameters produce 
approximately the same speed for engulfment. (B) Ranges of parameter values for successful engulfment. Simulation time was 
restricted to twice the engulfment time with Standard Parameters (SP) given in Methods and represented by vertial dashed 
line. Shown are fold changes of bending stiffness (blue), surface tension (green), ligand-receptor energy density (orange) and 
volume constraint (purple) along x-axis. Pointed arrows indicate that parameter regimes may extend beyond tested limits, 
blunt arrows indicate sharp limits beyond which full engulfment is not reached in simulation time. (C) Cross section of particle 
engulfment with particle radii R = 1.2/im (left) and R = 3.8/im (right). (D) Engulfment of a spheroidal particle with tip (left) 
and long side (right) flrst. Principal axis of the spheroidal particle are Ri = R2 = 1.5/im, and R3 = 4.2/im. Dashed curve 
indicates particle outline. (E) Stalled engulfment of spiral-shaped particle, characterised by a volume similar to a spherical 
particle of radius 2.2/im. 
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type energy function [55], including contributions from 
ligand-receptor binding, membrane bending and stretch- 
ing, as well as a cell- volume constraint. Chosen mem- 
brane parameters effectively describe the cell plasma 
membrane with its underlying actin cortex. The model 
was implemented using finite-temperature Monte Carlo 
simulations of the discretized cell membrane (see Meth- 
ods for details). Briefly, the algorithm proposes ran- 
dom, thermally generated membrane fluctuations (trial 
moves), which are either accepted or rejected depending 
on the change in total energy during the move. When a 
random membrane fluctuation brings the cell membrane 
in contact with the particle, this fluctuation is likely 
accepted due to energetically favorable ligand-receptor 
binding. Once accepted, this fluctuation is made irre- 
versible as a result of signaling and actin polymerization. 
In contrast, a membrane fluctuation far away from the 
particle is less likely to be accepted. Even if accepted, the 
fluctuation is not made irreversible and hence may retract 
at a later time (see Supplementary Fig. SI). Hence, in 
our model the actin network only supports membrane 
fluctuations which lead to progression of the engulflng 
zipper as a result of signaling. 



Dependence of phagocytic cup shape on membrane 
biophysical parameters 

Using our model for the zipper mechanism, we have 
successfully simulated phagocytic engulfment in a broad 
range of parameter values (see Figure 2). Figure 2 A 
shows two different characteristic cup shapes we ob- 
tained. Low surface tension (i.e. low energy cost for 
stretching the membrane and underlying actin cortex), 
and tight cell- volume constraint {i.e. high energy cost 
for increasing the cell volume) , lead to a thin phagocytic 
cup since a thin cup requires extra membrane but lit- 
tle extra volume. In contrast, weak volume constraint 
and/or high surface tension produce a broad cup. Based 
on the parameters explored, we chose intermediate val- 
ues for both surface tension and cell-volume constraint 
as our Standard Parameters (SP) for the remainder of 
our simulations in order to produce realistic cup shapes 
(see Methods for details). Figure 2B shows that most pa- 
rameter values can be changed independently by at least 
one order of magnitude, without negatively affecting en- 
gulfment completion. Note that changing simultaneously 
several parameters may affect engulfment more drasti- 
cally. Our simulations also show that cup shape depends 
on the kinetics of engulfment, determined by membrane 
fluctuations and therefore temperature (see Supplemen- 
trary Fig. S2). Additionally, preventing thermal fluctua- 
tions (by setting the temperature to zero Kelvin) during 
a simulation stops cup progression. This indicates that 
in our model membrane fluctuations are indeed required 
to bring receptors in close contact with ligand molecules 
on the particle, emphasizing their important role in the 



ratchet mechanism. 



Active versus passive uptake and the role of actin 

Although phagocytosis generally involves active pro- 
cesses such as actin polymerization in the cup (active 
zipper) [11], recent reports indicate that phagocytosis 
may still work in an act in-independent manner. In- 
deed, phagocytic uptake was observed despite treating 
phagocytes with (modest amounts of) cytochalasin D 
[39]. Hence, ligand-receptor binding may be sufficient 
in guiding the cell membrane around the particle under 
certain conditions (passive zipper). To investigate the 
energetic requirements of the zipper mechanism, we 
implemented simulations of the passive zipper. In these 
simulations, ligand-receptor bonds are not stabilized by 
actin polymerization and can unbind at later times, i.e. 
remain reversible. Hence, engulfment may still progress 
if the energetic cost of stretching and deforming the 
membrane is offset by the ligand-receptor binding energy 
in the presence of thermal membrane ffuct nations. 

Figure 3 {left) shows that engulfment of small (1.5/im 
radius) particles by the passive zipper leads to more 
variable phagocytic cup shapes than engulfment by the 
active zipper. For the active zipper, random membrane 
fluctuations are rectifled by irreversible ligand-receptor 
binding due to actin polymerization. This leads to 
uniform progression of the cell membrane all around 
the particle at approximately the same speed (Figure 
3A). In contrast, engulfment by the passive zipper 
occurs through binding of large membrane ruffles which 
eventually enclose the particle (Figure 3C). The vari- 
ability of the phagocytic cup may be a measure of the 
respective contributions of active and passive processes 
in engulfment progression. 

Figure 3 {right) compares time courses of the membrane 
energy and progression of uptake for the active and 
passive zippers. We found that in both cases, membrane 
energy decreases rapidly at the very beginning of the 
uptake process due to energetically favorable ligand- 
receptor binding without large, energetically unfavorable 
deformations of the cell membrane. After this short 
initial period, the total energy increases with simulation 
time. This increase is much more pronounced for the 
active zipper, which stabilizes energetically unfavorable 
random membrane fluctuations by actin polymerization. 
In contrast, for the passive zipper such high-energy 
membrane deformation may not last over time. The 
slower increase in energy for the passive zipper correlates 
with a slower engulfment. For the simulations shown, 
engulfment for the active zipper is approximately twice 
to three times as fast as for the passive zipper, although 
the latter eventually engulfs the particle. However, 
the difference between active and passive engulfment 
depends on biophysical parameters, and may be reduced 
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FIG. 3: (A) Side view (top) and cross section (bottom) of a phagocytic cup obtained for the active zipper. (B) Corresponding 
time courses of the membrane energy (thick solid, medium dashed and thin dotted black lines) and percentage of engulfment, 
defined by the average membrane height around the particle (blue solid, blue dashed, blue dotted lines) for three repeats of the 
stochastic simulation. (C) Side view (top) and cross section (bottom) of the phagocytic cup obtained for the passive zipper 
(same overall simulation time). (D) Corresponding time courses of the membrane energy (thick solid, medium dashed and thin 
dotted black lines) and engulfment (blue solid, blue dashed, blue dotted lines) for three repeats of the simulation. Dashed light 
blue line indicates corresponding maximal membrane height for one of the simulations (in percentage of the particle diameter). 
Small particles with 1.5/im radius were used. 



for particular choices of the parameters values, e.g. lower 
surface tension and/or stronger ligand-receptor binding. 



Particle size matters for passive, not for active 
zipper 

Experiments show that phagocytosis is relatively 
insensitive to particle size [42]. Using our model for 
the active zipper, we simulated engulfment of spherical 
particles with different radii ranging from 1.2 to 3.8/im. 
Figure 2C shows that engulfment progresses normally 
for small and large particles. Hence the active zipper 
mechanism is sufficiently robust to allow engulfment 
of differently sized particles using the same set of 
biophysical parameters, although engulfment of large 
particles requires more time. Noticeably, large particles 
(in general, with radius larger than 2.5yrxm) were taken 
up via more regular phagocytic cups than small particles 
(with radius 1.5/im or smaller), indicating that active 
processes may be required for engulfment of large 
particles. To confirm this observation we have simulated 



engulfment of large 3/im-radius particles by both the 
active and the passive zipper, shown in Supplementary 
Fig. S5. While the active zipper resulted in complete 
uptake of the particle, the passive zipper only engulfed 
a few percent of the particle's surface area. Thus the 
difference in phagocytic efficiency between the two zipper 
types was exacerbated for large particles, reflecting the 
importance of actin polymerization for engulfment of 
large particles. 



Experimental test of model predictions 

To test our model predictions and to specifically com- 
pare active with passive engulfment for small and large 
particles, we transfected COS-7 cells with either wild- 
type FC7R (WT-FC7R) or a signaling-dead mutant re- 
ceptor (Y282F/ Y298F-FC7R). Cells expressing the wild- 
type receptor are expected to perform active engulfment, 
whereas cells expressing the signaling-dead mutant re- 
ceptor are expected to perform passive engulfment. As a 
control, passive engulfment is additionally implemented 
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FIG. 4: (A) Typical fluorescence image of COS-7 cell with FC7R-GFP shown in green and IgG antibodies on 1.5/im-radius 
particles shown in blue. (B, C) Illustrative schematics of 3-dimensional confocal fluorescence microscopy images. (B) Schematic 
(side view) of a regular cup. (C) Schematic (side view) of a variable cup. From the intersections of the confocal imaging planes 
with the phagocytic cup, the maximal height of the cell membrane is determined for each angular segment around the particle 
(red dots). (D) Membrane height as a function of angular segment (see inset in panel F), corresponding to the cup drawn in 
panel B (dashed line) and C (solid line). Only the angular segments from index 1 to 12 out of 24 are shown. The variability 
of the cup is deflned as the standard deviation divided by the average of this function. The higher this measure, the more 
variable the cup. (E, F) Side and top views of half-engulfed particles, with cups reconstructed from confocal microscopy data. 
(E) Top view, taken in the equator plane of the particle, shows a regular distribution of cell membrane (receptors) around the 
particle (regular cup). (F) Top view shows an irregular distribution of cell membrane (variable cup). Inset Schematic deflnes 
the angular segments (in top view). 
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by treating cells expressing WT-FC7R with 0.2/iM of cy- 
tochalasin D (WT-Fc7R+CytoD), which prevents actin 
polymerization (see Methods). Synchronized phagocyto- 
sis assays using small (1.5 /im radius) and large (3 /im 
radius) IgG-opsonized polystyrene particles were carried 
out, and, after fixation, receptor localization in phago- 
cytic cups was visualized by fluorescence confocal mi- 
croscopy. Cells were imaged at different time points dur- 
ing phagocytosis, and at each time point, three to eight 
imaged cells were each engulflng simultaneously between 
four and twenty particles (Figure 4A). Consequently, for 
each condition we analyzed at least seventy particles. 

To test whether passive engulfment leads to more vari- 
able cups than active engulfment, we developed an image- 
analysis method illustrated in Figure 4B-F. The cup 
shape varibility was quantified by the standard deviation 
of the distribution of cell-membrane (FC7R-GFP fluores- 
cence) height around the particle, divided by the square 
root of the average membrane height. The unit of mem- 
brane height is given by the distance (0.4/im) between 
consecutive confocal image planes (see Methods). Figure 
5A shows that for small particles engulfed between 20 and 
40% of their surfaces, cup variability increases from cells 
transfected with wild-type receptor to cells transfected 
with signaling-dead mutant receptor to WT-Fc7+CytoD 
cells. The lowest variability, found for cells expressing 
wild- type receptor, is statistically signiflcant against both 
passive zipper types (Student's t-test, p- value < 0.001). 
This result is consistent with model predictions: Figure 
5 A, inset shows the cup variability from simulations, re- 
vealing that the active zipper leads to significantly less 
variable cups. In contrast, for the ranges of engulfment 
between 40 and 60% (Figure 5B) and between 60 and 
100% (see Supplementary Fig. S7) we observed no no- 
ticeable difference in cup variability between the three ex- 
perimental conditions, while our model consistently pre- 
dicts more variable cups for the passive zipper (Figure 
5B, inset). This discrepancy may indicate that active 
processes such as contraction by myosin motor proteins 
become important at later stages of engulfment, limiting 
our model's full validity to the early events in phagocy- 
tosis (see (7onc/^/szon section) . 

To confirm that our results are independent of the 
specifics of the analysis method used, we also analyzed 
phagocytic cups with an alternative, albeit less accu- 
rate, method (see Supplementary Fig. S9). This method 
consists of determining the distribution of FC7R-GFP 
fiuorescence intensity around the particle at its equa- 
tor plane, restricting the analysis to approximately half 
taken up particles (30 — 70%). The standard deviation of 
this distribution provides an alternative measure of the 
cup variability. We arrived at the same conclusion, con- 
firming our result of the difference in cup variability (see 
Supplementary Fig. 10). Finally, note that temperature- 
induced synchronization is imperfect and may lead to 
variability in cup growth [39]. However, our measures 



of the variability in cup shape are independent of such 
an effect since we include all time points together in the 
analysis (except for plots showing the time dependence 
of engulfment). 

Our experiments further show that cup shape ranges 
from regular to variable for all three experimental con- 
ditions, but that the frequency of different cup shapes 
depends on the condition. Figure 5C plots the reparti- 
tion of phagocytic cups for different conditions into both 
regular and variable cups. Note that a cup was identified 
as regular if its variability was below the variability aver- 
aged over all experimental conditions. In contrast, a cup 
was identified as variable if its variability was above the 
overall average. This plot shows that in our experiments 
a regular cup is most likely produced by a cell expressing 
wild-type receptor, whereas a variable cup is most likely 
to be produced by a cytochalasin-D treated cell. Hence, 
the cup shape has universal features independent of bio- 
chemical details. Examples of a regular and a variable 
cup are provided in Figures 4E and F, respectively. Both 
cups were taken from a cytochalasin-D treated cell, con- 
firming that a regular cup can occur under any of our 
experimental conditions. 

Our model also predicts that uptake by the active zipper 
is significantly faster than with the passive zipper (see 
Figures 3B and D). We experimentally tested this pre- 
diction by determining the percentage of engulfed surface 
area for each particle for different time points after ini- 
tiation of phagocytosis, and comparing this result with 
our simulations, in which simulation time was matched to 
actual time. Figure 6A shows that cells transfected with 
the wild-type receptor (active zipper) engulf significantly 
faster (three to four times) than cells under the other 
two conditions (passive zippers). This result is in quan- 
titative accordance with our model predictions (Figure 
6B). Furthermore, we determined the time dependence of 
phagocytic uptake for large particles. The active zipper, 
although slower for large than for small particles, still 
engulfs regularly, both in experiments (Figure 6C) and 
simulations (Figure 6D). Note that predicted and mea- 
sured time courses are in very good agreement without 
rescaling the time axis of the large-particle simulation. 
Furthermore, Figures 6C and 6D demonstrate the inabil- 
ity of the passive zipper to take up large particles, in 
both experiments and simulations. After more than 10 
minutes, the average engulfed surface area remains below 
20%. 

Note that for the time points measured, the average up- 
take does not exceed 50%, even for WT FC7R. This is 
caused by the fact that some particles are not engulfed 
irrespective of the condition [39] (see Supplementary Fig. 
S6), reducing the average percentage of engulfment. The 
proportion of almost completely engulfed particles (with 
engulfed surface area larger than 70%) beyond 6 min- 
utes is represented in the inset of Figure 6 A, showing 
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FIG. 5: (A) Experimental cup variability, defined in Figure 4, for 20-40% engulfed particles, for COS-7 cells transfected 
with wild-type FC7R (WT-FC7R) or signaling-dead mutant (Y282F/Y298F-Fc7R) in DMEM medium, as well as WT-FC7R 
cells treated with 0.2/iM cytochalasin D (WT-Fc7R+CytoD). The data points represent the variability averaged over all the 
phagocytic cups measured for given condition (receptor type and treatment) and range of engulfment. Error bars represent the 
statistical standard error of the cup variability. Student's t-test: p- values are (*) 0.0002 and (**) 0.0003. (Inset) Theoretical cup 
variability calculated from corresponding model simulations. Data points and error bars represent respectively the statistical 
average and standard error of the cup variability over several simulations. Student's t-test (*): p- value is given by 0.0002. (B) 
Experimental cup variability for 40-60% engulfed particles. (Inset) Theoretical cup variability calculated from corresponding 
model simulations. Student's t-test (*): p- value is given by 0.0001. (C) Repartitioning of regular and variable cups for different 
experimental conditions, i.e. wild- type receptor (WT-FC7R), signaling-dead mutant (Y282F/Y298F-Fc7R) and cytochalasin-D 
treated cells (WT-Fc7R+CytoD), for engulfment between 20 and 40%. The numbers indicate the contribution (in %) of the 
three conditions in the overall sample of regular (left) and variable (right) cups. 



that completion of engulfment is possible even for cells 
without act in polymerization. Note that long phago- 
cytic assays (12, 14, and 45 minutes) were performed 
for cytochalasin-D treated cells only, explaining why the 
difference in complete engulfment with cells expressing 
the wild-type receptor is smaller than one may expect. 

From these results, we conclude that functional act in 
polymerization is required for fast and regular engulf- 
ment in phagocytosis. Nevertheless, in line with our 



model predictions, cells showing deficient actin polymer- 
ization at cups are still able to uptake small particles, 
although more slowly and with more variable cups. 



Active zipper reproduces particle-shape dependence 
of phagocytosis 

Previously published experiments show that phagocy- 
tosis depends strongly on particle shape. In particular. 
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FIG. 6: Experimental data for small particles with 1.5/im radius (A) and large particles with 3/im radius (C). The data points 
represent averages of phagocytic cups obtained at given time points for the three different experimental conditions. Error bars 
represent the corresponding standard errors. (A, Inset) Proportion of almost completely taken up particles (engulfed surface 
area larger than 70%) after at least 6 minutes for the three conditions. (B, D) Progression of phagocytic engulfment with time 
in the simulations for small (B) and large (D) particles. Simulation time was converted to actual time so that 50% average 
engulfment of small particles by the active zipper corresponds to 6 min. No further adjustment of time conversion was done 
for large-particle simulations. The data points and error bars represent respectively the statistical average and standard error 
of phagocytic cups size from several simulations. 



elongated particles (similar to rod-shaped bacteria) are 
only taken up when presented to the phagocyte with their 
tip first [20, 38]. To test the particle-shape dependence 
of phagocytosis we conducted simulations of the active 
zipper using particles of different shapes while varying 
the initial orientation of the particle on the cell surface. 
Figure 2D shows the uptake of a prolate spheroid for 



two different orientations on the cell membrane after the 
same elapsed simulation time. In accordance with exper- 
imental observations [20, 38], uptake is more advanced 
for the spheroid particle engulfed with its tip first (about 
80% engulfed surface area) than for the particle attached 
along its major axis (about 50% engulfed in the same 
amount of time). This suggests a strong inhibitory effect 
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of high local curvature on uptake. In our model this is 
readily attributed to the energetic cost of bending the 
membrane around the two highly curved ends of an elon- 
gated particle placed horizontally on the cell membrane. 
In line with this explanation and experiments [38], the 
spiral-shaped particle in Figure 2E is not engulfed af- 
ter a time duration sufficient for the engulfment of large 
spherical particle of twice its volume. Hence, our simu- 
lations demonstrate that particle shape and orientation 
are indeed important biophysical parameters for phago- 
cytosis. 



Conclusion 

In this work, we studied the biophysical require- 
ments of the zipper mechanism, in particular the role of 
receptor-induced actin polymerization, and the effect of 
particle shape on uptake. In our model, the underlying 
biophysical mechanism of the zipper is an actin-driven 
thermal ratchet, which renders random membrane fluc- 
tuations irreversible close to the particle (Figure 1). This 
mechanism is supported by several recent experiments 
[37, 47, 54]. Previously, such Brownian ratchet models 
were successfully applied to explain force generation by 
actin polymerization and motility of the pathogen Lis- 
teria in hosts cells [49]. Our fully stochastic simulations 
can address for the first time the variability in particle 
uptake, recently noticed in phagocytic cup growth [39] 
and completion of endocytosis [46]. Implementation of 
our model in simulations led indeed to phagocytic engulf- 
ment for a broad range of values of membrane parameters 
(Figure 2), indicating exquisit robustness of the phago- 
cytic process. However, phagocytic cup shape depends 
on parameter values, specifically on the ratio between 
surface tension and cell- volume constraint, as well as on 
the kinetics of engulfment. Cells with low surface tension 
and/or tight volume constraint develop thin cups (Fig- 
ure 2 A, lejf), characteristic of Fc7R-mediated phagocy- 
tosis [4, 22]. In contrast, cells with high surface tension 
produce broad cups (Figure 2A, right). The latter cup 
shape is more reminiscent of CR3- mediated phagocytosis, 
although for this type of phagocytosis particles are be- 
lieved to sink into cells without protrusive cups [22, 23]. 

Using our model we were able to address the question 
whether the zipper mechanism requires an active driv- 
ing force, such as provided by actin polymerization. For 
this purpose, we compared the regular active zipper with 
a passive version of the zipper. In the passive zipper, 
ligand-receptor bonds are as strong as for the active zip- 
per (based on experimental observation [39]) but are re- 
versible, i.e. are not supported by actin polymerization. 
We demonstrated that the passive zipper also leads to en- 
gulfment of small particles (of radius 1.5/im), although 
cup progression is slower and more variable (see Figure 
3). Our active zipper can also reproduce the indepen- 
dence of uptake on particle size, in line with experimen- 



tal observations [42]. In contrast, large particles (of ra- 
dius 3/im) are poorly phagocytozed by the passive zipper. 
We subsequently confirmed these predictions with exper- 
iments by transfecting COS-7 fibroblasts with wild-type 
Fc7R and signahng-dead mutant Y282F/Y298F FC7R 
(see Figures 5 and 6). While the wild- type receptor 
represents the active zipper, the passive zipper is im- 
plemented through the use of signaling-dead mutant re- 
ceptor or treatment with cytochalasin D. Both prevent 
actin polymerization in the cups. Our study may indi- 
cate that ancient forms of phagocytosis were driven by 
physical (passive) principles, and only later in evolution 
biochemical regulatory pathways were added for further 
support and robustness. Passive phagocytosis may also 
become important when energy sources are scarce. 

Despite the robustness of phagocytosis to particle size, 
there appears to be a mechanical bottleneck around half- 
engulfment, recently observed by imaging [39] and also 
predicted by our model. For slight variations in some 
of the parameters, our simulations of the active zipper 
produce either complete or significantly incomplete up- 
take (see Supplementary Fig. S6). Indeed, when the cup 
grows, deforming the membrane costs more and more en- 
ergy per surface area engulfed. Beyond half-engulfment, 
the surface tension energy is twice as high as at the be- 
ginning of engulfment due to the membrane folding back 
onto itself. If this energy cannot be provided by the zip- 
per, then cup progression stalls and the particle remains 
incompletely engulfed. Alternatively, the experimental 
data on incomplete particle uptake may be the result of 
particle attachment on cell- membrane areas, unable to 
phagocytose due to other reasons, such as unfavorable 
local cell-surface curvature, proximity to cell edge or nu- 
cleus, or missing proteins, lipids, and smaller molecules 
belonging to key signaling pathways. Further studies will 
be required, ideally using live-cell imaging to avoid the 
need for conserving the cell's internal structures by fixa- 
tion [56]. 

Our model for the zipper mechanism can also explain 

the strong particle-shape dependence observed in phago- 
cytosis. Experiments show that an elongated spheroid is 
rapidly engulfed if the particle attaches to the cell mem- 
brane with its tip, but not if the particle attaches along its 
major axis [20, 38]. Furthermore, spiral-shaped particles 
are not phagocytozed [38]. The strong particle-shape de- 
pendence of phagocytosis is likely of biological relevance. 
On the one hand, it may increase the rate of infection of 
host cells by pathogenic bacteria. Indeed, recent experi- 
ments show that Mycobacteria tuberculosis and marinum 
are efficiently taken up and later released for spreading 
of the infection with the bacteria's tip first [57]. On the 
other hand, the highly curved shapes of some bacteria, 
e.g. the spiral-shaped Helicobacter and Campylobacter 
species, may prevent their uptake by macrophages [38], 
although injection of effector proteins can also be used 
by pathogens to hijack the immune or host cell's phago- 
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cytic response [58, 59]. Furthermore, the particle-shape 
dependence of phagocytosis may be exploited to improve 
drug delivery by enclosing active drugs in particles, whose 
shape prevents uptake and destruction by macrophages 
[38, 60]. 

While our biophysical model for the zipper mechanism 
is readily accessible for analysis and interpretation, the 
small number of model parameters makes a direct com- 
parison with measured parameter values difficult. First, 
our membrane parameters such as surface tension and 
bending stiffness are an order of magnitude smaller than 
reported bulk membrane parameters (see Methods). This 
reduction is not surprising as cells regulate these param- 
eters locally for efficient uptake [61-63]. Such regula- 
tion may include lowering of the surface tension by local 
membrane delivery through vesicles [35, 64, 65] and un- 
folding of membrane wrinkles [36], as well as changes in 
the lipid and protein composition in the phagocytic cup 
[1]. Second, our description of ligand-receptor interaction 
assumes that ligand and receptor distributions are con- 
tinuous and homogeneous, while experiments indicate the 
formation of receptor micro-clusters [66] , possibly as part 
of lipid rafts [47]. Third, our mechanism for actin poly- 
merization based on membrane fluctuations neglects the 
role of the motor proteins (myosin I and II), whose role 
in membrane deformations has been established [67, 68]. 
Consequently, our simulations describe well the dynamics 
of engulfment during the first two thirds of the uptake. 
At later stages of uptake, experiments show that phago- 
cytic cups close rapidly with a thin membrane protru- 
sion [4, 69], while our simulations show slow cup closure. 
However, taking into account myosin-driven contraction 
is beyond the scope of the current work. 

Our model may also be applicable to other biological 
systems in which a zipper-like mechanism is involved. 
One such example is sporulation of Bacillus subtilis 
during starvation. After asymmetric cell division, the 
larger mother cell engulfs the smaller forespore for spore 
maturation. Interestingly, the mother cell even engulfs 
the forespore when the cell wall is artificially removed. 
This process occurs in a fast, zipper- like fashion without 
known sources of energy [70]. Importantly, forespore en- 
gulfment is subject to high variation. About 60% of the 
cells successfully complete forespore engulfment, while 
40% do not at all, similar to the observation of the me- 
chanical bottleneck in phagocytosis. Other examples of 
engulfment may not be driven by a zipper mechanism. 
For instance, the penetration of red blood cells by the 
malaria Plasmodium merozoite is driven by an elaborate 
actin machinery of the particle (the parasite), which de- 
votes all its resources to wrap the cell membrane around 
itself. In stark contrast to phagocytosis, the engulfing 
host cell is completely passive [71]. 



Methods 

Theoretical techniques 

Cell membrane and particle models. The cell mem- 
brane is described by a two dimensional elastic sheet 
[55, 72], which includes both its lipid bilayer and as- 
sociated actin cortex. The particle is assumed to be 
rigid and immobile. Moreover, we account for ligand- 
receptor binding and include a constraint on the cell vol- 
ume. Therefore, the total free energy is given by 

E = + ^vol + ^LR, (1) 

where 

with (7^(r) = Ci^(r) + C2^(r) the square curvature ob- 
tained from the minimal (Ci) and maximal (C2) cur- 
vatures at point r. Note that the term corresponding 
to the product (7i(r)(72(r) is independent of the actual 
shape of the membrane as long as the overall topology is 
conserved, and therefore is ignored (Gauss-Bonnet the- 
orem). Bending stiffness hib refiects the energy cost of 
bending, and surface tension a reflects the energy cost of 
stretching the membrane with underlying actin cortex. 
Furthermore, expanding or shrinking the cytosol locally 
costs the energy 

Kol = ^cell(V) - ^cell(Vo) = ^p{V - Vq)', (3) 

where the quadratic dependence on actual volume V 
comes from the lowest order Taylor expansion of the 
cell energy around local steady state volume Vb, and 

K.p = ^ (^^f^) • Taylor expansion is justifled by 

the fact that our simulations and experiments use par- 
ticles significantly smaller than the cell (at least 10-20 
times in volume). Finally, in our model ligand-receptor 
binding is not described explicitly at the molecular scale, 
but is accounted for by a membrane-particle contact po- 
tential VLi?(i*), where Vlr{i^) = —Vlr if ^ membrane 
patch is within a distance Rq of the particle and zero if 
further away. Specifically, the associated energy is given 

by 

E^n= I dPvV^Kiv), (4) 

J m 

where the integral is performed over the cell-membrane 
area. Effectively, is given by the product of 

the individual ligand-receptor binding energy and the 
density of ligand-receptor bonds, divided by the density 
of vertices on the model membrane (see below). The 
width of the square potential Rq is chosen to be very 
small compared to the other length-scales of the model 
Rq < O.IR and does not infiuence the results. 
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Finite-element approximation. Simulations of 
phagocytic engulfment were implemented by discretizing 
the cell and particle surfaces using the Surface Evolver 
software [73] . This software is designed to perform energy 
minimization on flexible surfaces, and is freely available 
from http://www.susqu.edu/facstaff/ b/brakke/ evolver. 
The software includes a built-in programming language, 
which we used to implement a Monte Carlo algorithm 
(see below) . The cell membrane is approximated by a fi- 
nite number of vertices, used to create a triangular mesh. 
The software computes the local energy density at each 
vertex and sums up the energy contributions from all 
the surface elements to obtain the total free energy Eq. 1. 

Main model parameters. Our model uses four 
tunable biophysical parameters. Unless otherwise 
specified, we have used the set of Standard Parameters 
(SP), chosen according to experimental measurements 
[37, 42, 61, 74] when possible (see Supplementary Infor- 
mation, section 1), but ultimately to produce realistic 
cup shapes (see Figure 2). This set of parameters 
includes: the cell membrane bending rigidity and 
surface tension a, respectively set to 1.3 x 10~^ pN/xm 
and 6.2 x 10~^ mNm~^, i.e. slightly below the experi- 
mentally measured values since local changes in chemical 
composition of the cups membrane may reduce these 
parameters [1, 8, 36, 75-77]. The third parameter is the 
total binding energy density e = 58.5 pN/im~^. This 
value was estimated from measurements of the individual 
Fc7R-IgG binding free energy AFlr ^ 20kBT [78-81], 
the average density d^R = 270 — 435 /im~^ of IgG-Fc7R 
bonds [82], and the fact that in response to diffusion 
and trapping or signaling, receptors may cluster at the 
cup. Finally, the local constraint on cell-volume has 
been chosen Kp = 2.56 x 10~^ pN/im~^ to allow 20 
percent volume variation in line with observation [83, 84]. 

Monte Carlo simulations. The Surface Evolver was 
only used to obtain a triangular mesh (vertices connected 
by edges) of the cell membrane, and to resample the 
membrane as the uptake progresses. The cell- membrane 
evolution was implemented using finite-temperature 
Monte Carlo Metropolis simulations [85, 86]. Details 
of the simulation can be found in the Supplementary 
Information, section 2. Briefly, our algorithm calculates 
the total energy of the initial membrane configuration, 
then randomly selects a point to be the center of a 
membrane fluctuation, and a random direction and 
lateral extension of the fluctuation. Subsequently, the 
energy of the new membrane configuration is calculated, 
and compared to the initial energy. If the membrane 
fluctuation decreased the energy, the trial fluctuation is 
accepted and the procedure is reiterated starting from 
the new membrane configuration. To the contrary, if 
the membrane energy increased with the trial fiuctua- 
tion, the latter may be rejected with some probability 
depending on the configuration's energy difference. 
In this case, a new fluctuation is attempted from the 



initial conflguration. Between trial fluctuations, the 
cell-membrane vertices are examined. For the active 
zipper, the vertices within the closed neighborhood of 
the particle are immobilized for the remainder of the 
simulation. For the passive zipper, every membrane 
fluctuation may be reversed at a later time. 

Additional simulation parameters. (1) The ampli- 
tude of a membrane fluctuation. This parameter is set 
to 0.5i?o in all the simulations. (2) The mesh size, de- 
scribing the maximal distance between two neighboring 
vertices. This parameter is set to Rq. (3) Minimal width 
of a membrane fluctuation. To ensure that a minimum 
number of vertices is involved in each fluctuation, this 
parameter is set to 4i?o- (4) Mesh reflnement range. 
Our simulation script automatically refines the mesh lo- 
cally around a previously immobilized vertex within this 
range. This parameter is set to the particle radius i?, 
leading to reasonably smooth cup shapes in a reason- 
ably short calculation time (1-2 days for complete up- 
take using a Intel(R) Core(TM)2 Quad CPU working at 
2.50GHz and run by the RedHat EL5 linux distribution 
with 4GB RAM). 



Experimental techniques 

Cells, plasmids and antibodies. COS-7 fibroblast 
cells were obtained from American Type Culture Col- 
lection (ATCC) and cultured in Dulbecco's Modified 
Eagle's Medium (DMEM), supplemented with 10% foetal 
bovine serum (FBS) and penicillin/streptomycin (Invit- 
rogen). cDNAs encoding wild-type and Y282F/Y298F 
human Fc7RIIa from pRK5-Fc7RIIa [87] and pRK5- 
Y282F/Y298F-Fc7RIIa [88] were subcloned into pEGFP- 
Nl (Clontech) using primers 5'-ggtccaactgcacctcggt-3' 
and 5'-ccccccgaattctgttattactgttgacatggtc-3'. The 
cytoplasmic tail truncation mutant 239-Fc7RIIa 
was generated from the pRK5-Fc7RIIa template 
using primers 5'-ggtccaactgcacctcggt-3' and 5'- 
gggggggaattctcctgcagtagatcaaggccact-3 ' . Rabbit 
anti-bovine serum albumin (BSA) serum was purchased 
from Sigma- Aldrich. Alexa-conjugated secondary anti- 
bodies and phalloidin were purchased from Invitrogen. 

Transfection and phagocytic challenge. COS-7 
cells were transfected with GFP-tagged Fc7RIIa con- 
structs using an Amaxa Nucleofector and Nucleofector 
cell line kit R following the manufacturer instructions. 
For phagocytosis assays, transfected cells were seeded 
onto glass coverslips in 24-well plates at a density of 
15,000 cells/coverslip and incubated at 37 °C for 72 h. 
1 hour before commencement of phagocytosis assays, 
cells were incubated for 1 hour at 37 °C with serum- free 
DMEM plus 10 mM Hepes (Invitrogen). 1.5/im- and 
3/im-radius latex-polystyrene particles (Sigma- Aldrich) 
were opsonized by first incubating overnight at 4°C with 
3% BSA fraction V in PBS (Sigma- Aldrich) followed 
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by incubation with 1:100 dilution of rabbit anti-BSA in 
PBS for 1 hour at room temperature. Particles were 
re-suspended in ice-cold serum-free DMEM plus 10 mM 
Hepes at a concentration of 1.5x10^ particles/ml and 
500 /il added to each coverslip. Plates were incubated 
on ice for 10 minutes to allow binding of particles. 
Medium was then replaced with pre-warmed serum-free 
DMEM plus 10 mM Hepes and plates were incubated 
at 37 °C, then processed for scoring or microscopy 
as described below. Experiments carried out on cells 
treated with cytochalasin D were conducted as above 
with the addition of a 20 min pre-incubation step with 
0.2 jjM cytochalasin D in serum-free DMEM plus 10 
mM Hepes at 37 °C immediately before incubation of 
cells with opsonized particles. This concentration of cy- 
tochalasin D was included in all further incubation steps. 

Imaging of phagocytic cups. Plates were placed on 
ice after a 20 minute incubation at 37 °C and medium 
replaced with a 1:500 dilution of anti-rabbit Alexa 488 
in 3% BSA/PBS at 4°C for 5 min. Cells were fixed after 
incubation at 37 °C for the appropriate amount of time 
with ice-cold 4% paraformaldehyde/PBS, permeabilised 
and labelled with goat anti-rabbit Alexa 633, and 
phalloidin Alexa 555 for visualizing F-actin at room 
temperature for 30 min. Z-series image stacks were 
acquired on a Zeiss LSM-510 confocal microscope using 
a step size of 0.4 /im. 

Image analysis of phagocytic cups. Two fluorescence 
channels (IgG, FC7R-GFP) were acquired and analyzed 
using MATLAB (MathWorks). During the acquisition 
process, the FC7R-GFP fluorescence intensity was set to 
zero at any pixel where the IgG intensity is null. The 
fluorescence intensity distribution of IgG was used to 
determine the coordinates of the centers of particles with 
their corresponding radii, using an automated search 
based on the Hough transform [89], available online at 
www.mathworks.com/matlabcentral/flleexchange/4985 
within each 2D image. The percentage of engulfed par- 
ticle surface area was calculated by comparing the local 
FC7R-GFP (cell membrane) and IgG (particle surface) 
intensity distributions within a sphere Sq of radius AR/S 
whose center coincides with the particle center. Two 
methods were used to quantify the variability of the 
cups. The method used in the main text cuts the three- 
dimensional image (and hence the circular projections 
of particles within each imaging plane) in twenty-four 
angular segments, and flnds the highest plane in which 
FC7R-GFP fluorescence intensity is detected in the 
immediate neighborhood of the particle for each segment 
(see Figure 4). This analysis produces a distribution 
of membrane height versus angular segment index, for 
which we compute the average and the mean-square 
deviation. The average height reached by cell mem- 
brane is roughly proportional to the surface engulfed 



(see Supplementary Fig. SIO), and the mean-square 
deviation divided by the square root of the average 
height quantifles the variability of the phagocytic cup, 
excluding a trivial size dependence on the variability. 
The less accurate alternative method determines the 
angular distribution of FC7R-GFP fluorescence intensity 
within the particle's equator plane in the particle's 
immediate neighborhood, keeping only particles whose 
uptake level comprises between 30 and 70% (roughly 
half-engulfed particles). Speciflcally, this region is cut 
into twenty- four identical angular segments, for which 
the total FC7R fluorescence intensity is calculated. Then 
the average intensity per segment is calculated, as well 
as the standard deviation. The higher the standard 
deviation, the more variable the cup. 



Accession Numbers 

The Fc7RIIa receptor is referenced in protein database 
Genbank (http://www.ncbi.nlm.nih.gov/Genbank) un- 
der the accession number CAA01563. Signaling-dead 
mutant receptor Y282F/Y298F-Fc7RIIa is obtained by 
replacing tyrosines (Y) with phenylalanines (F) at posi- 
tions 282 and 298. 
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Additional discussion on parameter values and 
simulations 

1 - Choosing model parameters 

According to previous experimental studies of phagocytes, 
a cell's bulk bending rigidity Kb may vary in a range from 
0.03 — 2 pN/^m depending on the type of cell and measurement 
method used (see [1] and references therein). Similarly, the 
bulk surface tension a may range from 2 — 4xl0~^ mNm~^ 
in resting neutrophils but can vary in other cell types [1- 
3]. However, during phagocytosis of large particles with ra- 
dius 5/im by aspirated spherical neutrophils, an abrupt in- 
crease in surface tension was observed (to 0.8 mNm""*^ [3]) 
during the uptake. Note also that, hit and a may vary on 
the scale of the cell [4-6]. Changes in local phospholipid 
and protein composition of the cups membrane may influence 
Kb [7]. Furthermore, during engulfment membrane delivery 
from internal buffers [8-11] and unfolding of surface mem- 
brane folds [12] may strongly reduce a at the cup. Therefore 
in our simulations we use values of Kb and a slightly lower 
than previously cited values, i.e. Kb = 1.3 x 10~^ pN/xm 
and a = 6.2 x 10~^ mNm~^ in order to produce realistic cup 
shapes. 

The Fc7R-IgG binding free energy AFlr was measured to 
WksT [13], 18-20/cbT [14, 15], and 23k bT for IgE to Fee RI 
binding [16] at physiological temperature. The density cIlr of 
IgG-Fc7R bonds was estimated to range from 270 — 435 fim~^ 
for macrophages [17]. The total binding energy density e is 
consequently of the order of AFLndLR = 30 pN/im~^. Note 
that due to simple diffusion or in response to signaling, addi- 
tional receptors may be recruited to the phagocytic cup, which 
may lead to a higher binding energy density. Additionally, re- 
ceptors may cluster, leading to a non-uniform ligand-receptor 
binding-energy density in cells. Taken together, we chose 
a ligand-receptor binding energy density e = 58.5 pN/im~^ 
slightly higher than measured in resting cells. 

The cell- volume constraint was chosen Kp = 2.56 x 
10~^ pN/im~^. This reflects strong regulation of cell volume 
on the one hand [18, 19], but, on the other hand, allows vol- 
ume fluctuations of about 20%, in line with observed distri- 
butions of cell volumes [20, 21]. 

Finally, the algorithm convergence is maximized if approx- 
imately half of the trial moves are accepted. To satisfy this 
requirement without constraining the biophysical parameters, 
we have modulated the trial moves acceptance rate by choos- 
ing an intermediate cutoff i^max for the width of the mem- 
brane fluctuations (see below). The above listed parameter 
values for Kb, cr, e and Kp constitute the standard set of pa- 
rameters values. 



2 - Details of Monte Carlo algorithm 

Within the Surface Evolver, we implemented finite- 
temperature simulations using a Monte Carlo Metropolis al- 
gorithm [22, 23] to describe the membrane dynamics. Our 
Monte Carlo algorithm includes four steps: 
(1) Calculation of the total energy of the initial membrane 
configuration Ei. (2) Random selection of a vertex for a 
membrane fluctuation (trial move) and a vector u for the di- 
rection of the trial move. The lateral width w of the trial 
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FIG. 1: Ratchet model in action. Panels (A) and (B) show 
phagocytic uptake at two different times with Ia < tp for 
the same simulation. (A) The fluctuation (1) occurs close to 
the particle, whereas the fluctuation (2) does not bring the 
cell membrane in contact with the particle. Panel (B) shows 
that the membrane at position (1) has progressed around the 
particle, whereas the membrane at position (2) has almost 
completely retracted to its initial fiat state, demonstrating 
the importance of membrane immobilization by effective actin 
polymerization during uptake. 



move is randomly chosen between and an upper limit ic'max, 
possibly determined by the spacing between actin cortex- 
membrane linkage proteins [24]. The trial move consists in 
moving the vertex and all the vertices within a distance 
d < w from by the vector 0.5 (1 + cos{d/w)) u. This par- 
ticular function was chosen because of its regular properties 
at d = w (vanishing of the function and its first derivative). 
(3) Calculation of the total energy E2 after the trial move. (4) 
Decision to whether accept or reject the trial move is based 
on the Metropolis criterion. If E2 < Ei, the move is accepted 
because the new membrane configuration has a lower energy 
than the initial one. If E2 > Ei, the move is accepted with a 
probability exp [—{E2 — Ei) / (kpT)]^ where T is the simula- 
tion temperature in Kelvin. If the trial move is rejected, a new 
trial move is selected from the same initial configuration; if it 
is accepted a new trial move is selected from the new mem- 
brane configuration. Between each trial move, the cell mem- 
brane vertices are examined. For the active zipper, the ver- 
tices within a cutoff distance from the particle (chosen 0.5i?o) 
are immobilized, i.e. these vertices do not move anymore. As 
a consequence, the corresponding patch of membrane is sta- 
bilized, and the ligand-receptor bonds are irreversible. 

For the passive zipper, none of the vertices are immobi- 
lized and all the trial moves may be reversed at a later time. 
Note that our choices of simulation parameters imply that the 
number of accepted and rejected fluctuations are of the same 
order of magnitude, which improves the convergence of the 
Metropolis algorithm. 

Our equilibrium Monte Carlo approach is justified a pos- 
teriori to study the dynamics of uptake, since only a very 
small proportion (less than 1%) of the accepted trial moves is 
irreversibly immobilized. Hence, the algorithm samples mem- 
brane configurations extensively in line with equilibrium ther- 
modynamics. 
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FIG. 2: Relaxation of phagocytic cup to the ground state. (A) Phagocytic cup shape at half-engulfment obtained for the 
active zipper with the set of Standard Parameters (left); cup relaxed at zero absolute temperature (right). (B) Phagocytic cup 
shape at half-engulfment obtained for the passive zipper with standard set of parameters (left); cup relaxed at zero absolute 
temperature (right), leading to unzipping. 




FIG. 3: Influence of the maximal width of membrane fluctua- 
tions on cup shape. Shown are (top) narrow fluctuations with Q j 
maximal width Wma^ = 1.5i?, (middle) medium wide fluctu- % 
ations with Wma^ = 3i?, and (bottom) wide fluctuations with i 
i^max = 4i?. Parameter R is the particle radius. f 



Additional simulations 

3 - Ratchet model in action 

Our ratchet-like biophysical model is based on effective 
actin polymerization to render ligand-receptor bonds irre- 
versible. As ligand-receptor binding and hence signaling can 
only occur in the immediate neighborhood of the particle, 
random membrane fluctuations far from the particle may re- 
tract with increasing simulation time. Figure SI gives an 
example of a reinforced membrane fluctuation near the parti- 
cle and a vanishing membrane fluctuation further away from 
the particle, illustrating how this concept is implemented in 
the simulations. 



4 - Role of engulfment kinetics in cup shape 

Our simulations show that cup shape depends on the ki- 
netics of phagocytic uptake. To demonstrate this we re- 
laxed a phagocytic cup obtained with a flnite-temperature 
Monte Carlo simulation towards the ground state, i.e. the 



FIG. 4: Reproductibility of cup shapes for simulations under 
the same conditions. Side views of three simulations with the 
set of Standard Parameters (for approximately 50% engulf- 
ment) . 



shape at zero absolute temperature in the absence of ther- 
mal fluctuations (see Figure S2). During relaxation, the 
membrane evolves towards the minimal energy conflgura- 
tion (right) starting from the initial conflguration (left). In 
the case of the active zipper (panel A), the resulting cup 
not only smoothes, but also changes shape (here it becomes 
broader). For other parameter choices, simulated cups are ei- 
ther broader or thinner than predicted by their ground state 
due to flnite temperature stochastic engulfment, emphasizing 
the importance of the randomness of membrane fluctuations 
in shaping the cup. 

In the case of the passive zipper (panel B), flnite tempera- 
ture stochastic processes are even more crucial for engulfment. 
Indeed, Figure S2B (left) is obtained with a simulation of 
the passive zipper at flnite temperature. After some time, the 
simulation is stopped and restarted at zero absolute tempera- 



4 




Simulation steps (x 1000) 
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FIG. 5: Active versus passive engulfment of a large particle. (A) Cross section of phagocytic cup obtained for the active zipper. 
(B) Time course of total energy (black) and percentage of engulfment (blue). (C) Cross section of phagocytic cup obtained 
for the passive zipper (same overall simulation time). (D) Time course of total energy (black) and engulfment (blue) for the 
passive zipper. The particle radius is R = 3/im. Three simulations were carried out in (B,D). 




FIG. 6: Stalling at partial engulfment. (A) Time course of the percentage of engulfment as a function of increasing surface 
tension in simulations. Shown are completed (or nearly completed) uptake (black lines) and partial (stalled or very slowly 
progressing) uptake (blue lines). Surface tension values (cr x 10~^ mNm""*^) used are 0.8 (black solid), 1.1 (black dashed), 1.4 
(thick blue dash-double dotted), 1.6 (thick blue dotted), 1.9 (thick blue dashed-dotted), 2.1 (thick blue dashed), 2.3 (thick blue 
solid), 2.8 (thin blue dotted), 3.1 (thin blue dash-dotted), 3.4 (thin blue dashed), and 3.7 (thin blue solid). (B) Distribution 
of engulfed particle surface area measured in our experiments for COS-7 cells transfected with wild-type receptors (WT-FC7R) 
after lOmin. Black and blue bars correspond to the two peaks of the bimodal distribution. 
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FcyR: WT Y282F/Y298F WT 
Treat.: DMEM CytoD 

FIG. 7: Statistical analysis of phagocytic cup shapes. Cup 
variability (for 60 to 100% engulfed particles) was determined 
for COS-7 cells transfected with wild-type FC7R (WT-FC7R) 
or signahng-dead mutant (Y282F/Y298F-Fc7R), as weU as 
WT-FC7R cells treated with 0.2/iM cytochalasin D (WT- 
Fc7R+CytoD). The data points represent the variability av- 
eraged over all the phagocytic cups measured for given condi- 
tion and range of uptake. Error bars represent the statistical 
standard error of the cup variability. (Insets) Theoretical cup 
variability calculated from corresponding model simulations. 
Student's t-test (*): p- value is less than 0.0001. 




Surface engulfed (%) 



FIG. 8: Average cup height < h> versus percentage of en- 
gulfed particle surface area S. Each circle corresponds to an 
imaged cup (black for WT-FC7R, green for Y282F/Y298F- 
FC7R, and red for WT-Fc7R+CytoD). For a perfectly regu- 
lar cup, both quantities are related by 5'(%) = 100 < h > 
(/im)/(7 X 0.4). Due to shape variability and finite number of 
angular segments used, values obtained lie slightly above the 
diagonal straight line. 



area. The fiuctuation-acceptance rate is decreased, the cup 
grows more slowly and becomes broad (bottom panel). 



ture for relaxation of the partially engulfed state towards the 
ground state. After partial relaxation, the cup presented in 
Figure S2B (right) is obtained. The membrane smoothes 
similar to the active zipper, but also retracts, resulting in 
unzipping. This indicates that the release of ligand-receptor 
binding energy is not sufficient for uptake by the passive zip- 
per. Additional thermal membrane ffuctuations are also re- 
quired. Further relaxation eventually leads to a total retrac- 
tion of the cell membrane for the passive zipper. 



6 - Reproductibility of cup shapes 

Phagocytic cup shape is variable, in particular without 
the support of the actin cytoskeleton for the passive zipper. 
However, even for the active zipper, engulfment is inherently 
stochastic due to random membrane ffuctuations. As depicted 
in Figure S4, simulations with the same set of parameters 
lead to slightly different cups after the same simulation time. 
This illustrates particle-to-particle variation. 



5 - Influence of membrane fluctuation width on cup 
shape 

The maximal width of the ffuctuations i^max is chosen to 
optimize the convergence of the Monte-Carlo Metropolis al- 
gorithm (see Methods in the main text). This algorithm accu- 
rately samples the probability distribution of the membrane 
conffguration, and consequently of cup shapes at different 
stages of uptake. The convergence of the algorithm is op- 
timal if about half of the randomly selected membrane ffuc- 
tuations (trial moves) are accepted. This condition is fulfflled 
for intermediate values of lOmax- Figure S3 shows how the 
choice of Wma^ inffuences the cup shape. If Wma^ is small, only 
laterally small membrane ffuctuations are generated by the al- 
gorithm. These ffuctuations do not change the surface area 
and cell volume signiffcantly, and consequently are frequently 
accepted by the algorithm. Hence, the cup grows fast and 
thin (top panel). In contrast, if lOmax is large, laterally wide 
ffuctuations may be selected, and are less frequently accepted 
since they can lead to drastic changes is cell volume and/or 



7 - Active versus passive engulfment of large 
particles 

In the main text we demonstrate that the passive zipper 
is able to take up small particles of radius 1.5/im, although 
slower and in a more variable fashion than the active zipper. 
Here, we compare the ability of the two zippers to take up a 
larger particle of radius 3/im. Figure S5 shows that the pas- 
sive zipper only engulfs a few percent of the particle surface 
area (panels C and D), whereas the active zipper completes 
uptake (panels A and B). This demonstrates that, for the 
uptake of large particles, ligand-receptor binding in conjunc- 
tion with thermal ffuctuations is not sufficient, and that the 
functionality of the actin network is fully required. 



8 - Stalling at partial engulfment 

A model previously published by one of the authors uses the 
force generated by actin on the cell membrane as a parameter. 
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Angular segment Average fluorescence Fq 



FIG. 9: Alternative image analysis of phagocytic cups. (A) Side view of a typical 3-dimensional fluorescence image, recon- 
structed from confocal microscopy data. Shown are the FC7-GFP receptors (green) and IgG ligand coating the particle (blue). 
The analysis plane (white square) intersects the particle at the equator. (B) Top view of the analysis plane. The green light 
distribution is determined using the total green light intensity in the different angular segments. (C) An example of such a 
distribution for one particular particle, used for extracting the average intensity Iq and the mean-square deviation 51. (D) Illus- 
trative scatter plot and classification of particles with respect to their engulfment statistics Iq (horizontal axis) and normalized 
root mean-square deviation 61^ = 5I/\/Iq (vertical axis). Regular cups correspond to a smaller root mean-square deviation of 
the fluorescence distribution than variable cups. 
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FIG. 10: Alternative statistical analysis of phagocytic 
cup shapes. Cup alternative variability (arbitrary units) 
for 30 to 70% engulfed particles, for conditions WT- 
FC7R, Y282F/Y298F-FC7R, and WT-Fc7R+CytoD. The data 
points represent the variability averaged over all the phago- 
cytic cups measured for given condition and range of uptake. 
Error bars represent the statistical standard error of the cup 
variability. Student's t-test: p- values is 0.003 (*). (Insets) 
Whisker plot representation of the cup-shape alternative vari- 
ability. 



and predicts that slightly different force strengths lead to very 



different behaviors of the phagocytic cup [25]. Specifically, a 
bimodal distribution of phagocytic cups was found, i.e. cups 
either fully take up the particle or stall before 50% engulf- 
ment. This was interpreted as a mechanical bottleneck at half- 
engulfment (see main text). Our simulations also exhibit this 
effect. Figure S6A shows how slight (much less than an or- 
der of magnitude) variations in biophysical parameters (here 
surface tension) may have drastic effects on the completion 
of uptake. This high sensitivity to small parameter changes, 
for specific parameter ranges, explains the particle-to-particle 
variation observed in our experiments (Figure S6B). 



Additional analysis of fluorescence images 

9 - Statistics of phagocytic cup shapes 

In Figure 5 of the main text, we demonstrate that the cup 
variability in experiments is increased when signaling to the 
actin cytoskeleton is interrupted or when cells are treated with 
cytochalasin D. More precisely, this result holds at the begin- 
ning of uptake between 20 and 40% of engulfed particle area. 
Between 40 and 60% of the uptake, there is no significant dif- 
ference between the three experimental conditions. For com- 
pleteness, we report here the cup variability between 60 and 
100% of uptake. Figure S7 shows that almost completed 
cups are more regular (the variability is reduced by a factor 
of 2 compared to the range 40 to 60%), and we observe no 
statistically significant difference between the three different 
conditions. Note that the variability obtained in simulations 
(see Figure S7 inset) is also reduced by a factor of 2, and 
that the relative difference in variability between the active 
and passive zipper is also reduced, as expected since fully 
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completed cups are necessarily regular. 



10 - Average cup height versus percentage of 
engulfed particle surface area 

Our image analysis method makes use of two different mea- 
sures of uptake. The height h reached by the cell membrane 
(measured by confocal slice index), and the percentage of par- 
ticle surface S covered by cell membrane. The height is used 
to calculate the cup variability, and the surface engulfed is 
used to classify the cups (20 to 40%, 40 to 60%, and 60 
to 100%). Figure S8 shows that these two quantities are 
indeed highly correlated. For a perfectly regular cup, the 
height < h >, averaged over angular segments, is given by 
< h >= 75/100 X 0.4/im, because each particle is imaged in 
7 confocal planes. As shown, this relation does not exactly 
hold for variable cups, and the average height is slightly higher 
than predicted by the engulfed surface area. This minor de- 
viation is due to the finite number of angular segments used. 



11 - Alternative image analysis of phagocytic cup 
shapes 

To rule out any bias in the analysis, we additionally use 
an alternative method to characterize cup variability. The 



method is based on the observation that a perfectly regular 
cup (at half-engulfment) reaches the equator plane all around 
the particle circumference, corresponding to an engulfed sur- 
face area S = 50%, and an average height <h>= 1.5/im. For 
such a cup the distribution of FC7R-GFP fluorescence signal 
(indicating the presence of the cell membrane) should be uni- 
form around the projection of the particle onto the equator 
plane (circle of radius R). In contrast, a variable half-engulfed 
cup has an irregular distribution of fluorescence intensity as 
the particle is only partially covered by cell membrane in its 
equator plane. To exploit this observation, we calculate the 
distribution of FC7R-GFP fluorescence along the particle's 
perimeter in angular segments. The standard deviation of 
this distribution quantifles the alternative measure of the cup 
variability, illustrated in Figure S9. 

Figure SIO shows the result from the statistical analy- 
sis of phagocytic cup shapes as obtained with the alternative 
method. Although this analysis method is more restrictive 
and inaccurate, it conflrms our result from the main text. 
Using 30 to 70% engulfed particles we obtain that cells trans- 
fected with wild-type FC7R produce signiflcantly less variable 
cups than cytochalasin D treated cells (Student's t-test: p- 
value= 0.003). 
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